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^ollowirig  is  a  translation  of  an  article  by  V.  I. 

Gubin  and  S,  Karimberdyyeva  in  the  Russian-lan^age 
periodical  Izvestiya  Akadeinii  nauk  UzSSR»  Seriye. 
fiziho-niatematicheskiKh  nauk  (Bulletin  of  the  Acad- 
en^f  Sciences  Uz^k'  SSR“  PhysicomtheKatical 
Science  Series),  Tashkent,  No.  3^  19^0,  pages  38- 

^3.7 

In  connection  with  the  present  development  of  computer  tech¬ 
nology  there  are  great  potentialities  for  numerical  prediction  of 
the  weather  and,  in  particular,  for  the  F^-ediction  of  the  synoptic 
situation.  We  shall  examine  one  of  the  systems  of  the  precalcula¬ 
tion  of  the  geopotential  hei^t  of  the  700-mb  isobaric  surface  by 
means  of  the  "Ural”  electric  computer,  which  is  being  used  in  the 
calculation  center  of  the  Institute  of  Mathematics  of  the  Academy 
of  Sciences  Uzbek  SSR, 


Basic  Ratios 

For  the  calculation  of  the  geopotential  height  of  the  700^mb 
isobaric  surface  we  used  a  barotropic  model  of  the  atmosphere, 
though  in  this  system  a  change  of  the  height  of  the  isobaric  s^face 
is  determined  only  by  the  transfer  of  the  vortex,  it  is  nevertheless 
interesting  to  explain  the  effect  of  the  steric  and  time  spacings 
upon  the  quality  of  the  prediction,  namely,  on  this  simplest  model. 
The  original  equatioh  in  the  one-layer  problem  is 

A0  «  -  ^  +1)# 


where  §  is  the  geopotential, 

1  is  the  Kariolis  parameter, 
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Inasmuch  as  ^  =  gH  (g  =  9*8  m/sec^,  H  being  the  height  of  -Uie  iso- 
baric  surface),  then  equation  (l)  may  h'e  Rewritten  in  the  following 
form 


Ah  «  -  (H,  g  ^  +  1). 

•L 


We  shall  solve  this  equation  by  the  terminal  differences  method. 
For  the  solution  ve  chose  a  grid  field  with  a  spacing  Ss  =  25O  km, 
the  number  of  points  being  equal  to  39I  (where  x  is  23  points  and 
y,  17),  Equation  (2)  is  solved  in  the  condition  of  the  rotation 
of  H  to  zero  on  the  boundary  of  the  grid  field.  Let  us  introduce 
the  sjonbol 


e  «  (H,  g 


(3) 


Then  for  the  arbitrary  unit  with  index  i,  j  may  be  vrritten 


©Ij  =  I  (H,  AH)ij  +  (H,  l)i_3 
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Using  centered  differences,  we  shall  obtain  approximate  values  of 
the  derivatives : 
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The  expression  for  6,  if  (5)  is  taken  into  consideration,  will 
take  the  following  form; 

^1  IT 

®ij  ”  -  ^1,0-2  * 

Hi-l,j+l  -  •  (Hi+i^j  -  Hi-i^j)  - 
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Here 


g 

^  ^2  “  2  , 

.  ■  .• 

If  there  is  taken 


1-4  1  . 


1  «  1.2  •  10-^^  Juj  Iv  “  lo92 

sec  ^ 


10' 


1  ;  g  *  9.8 

meter  second  sec 


Ss  »  2*5  •  105  meter,  then  k;^  ^2  'will  have  the  values 

k-.  =  3.92  •  10-6 - 1 - kp  =  2  •  10-8  _1_^ 

sec  •  gpdkm  sec 


Let  us  Introduce  the  symbol 
®lj  “  ®ijS  s2. 


Then  from  (2),  representing  the  left  side  in  a  different  form,  we 
shall  obtain  for 
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Erom  the  ratio  (7)  it  is  possible  to  find  the  values  for  H  in  all 
units,  using  the  iteration  process.  However,  the  problem  may  be 
solved  more  quickly  (in  the  sense  of  the  use  of  mchine  time)  with 
out  a  decrease  of  accuracy,  if  a  local  system  is  used.  For  this 
we  will  consider  that  the  effect  of  the  second  series  of  points 
(with  respect  to  i,  j)  is  insignificantly  small,  and  the  contribu¬ 
tion  of  these  sources  of  0‘  can  be  disregarded.  Such  an  assump¬ 
tion  is  adequately  substantiated  by  the  presence  of  a  more  or  less 
uniform  distribution  of  positive  and  negative  sources  on  this 
local  form. 

In  our  example  of  caieulation,  nine  points  were  chosen  as 
units,  in  which  the  effect  of  the  sources  was  taken  into  considers 
tion:  point  i,  j,  in  which  is  being  sought  Hij,  and  eight  points, 
which  svirround  it,  -i+l,j;  i,j+l;  i-l,oj  i+l,j+l;  i-^j+1; 

i-l,u-lj  i+l,j-l.  Then,  setting  up  the  nine  equations  ly, 
taking  into  consideration  the  conditions  in  the  local  contour,  we 
shall  write  Hij  in  the  following  form; 
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We  dlso  used  this  expression  for  the  determination  of  a  change  of 
the  height  of  the  isobarlc  surface  in  each  unit  of  the  grid  field. 
Selefcting  the  spacing  by  time  At,  it  is  possible  to  then  deter¬ 
mine  the  next  field,  using  the  Euler  method. 


♦  H*J  .  At. 


(9) 


In  the  calculation  of  the  height  of  the  700-mb  surface  we 
used  At,  equal  to  1  hour,  1»5  hours,  and  2  hours o  After  the  fiei^ 
through  At  was  found,  the  values  of  new  H  were  substituted  in  (6) 
and  the  whole  ca3-culation  procedure  was  repeated.  Inasmuch  as  the 
prediction  of  the  field  was  given  for  2h  hours,  the  number  of 
spacings  by  time  was  equal  to  24,  l6,  and  12, 

In  the  calculation  process  errors  are  accumulated  which  can 
lead  (as  in  any  kind  of  differential  system)  to  the  appearance  of 
spupioiis  ’M'£LV0S  f  III  oi'd.^r  *fco  avoid.  *t]i0SG  spuxious  waves  ^  we  car— 
ried  out  a  smoothing  of  the  field  every  6  hours  according  to  for- 
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Hsm 


I  Sij  +  (Hi+1,  j  +  %,d+l  +  Hi-l,d  +  Hi,  j-1^ 


(10) 


The  ratios  (6),  (8),  (Q),  and  (lO)  are  basic  for  the  solution  of 
the  problem  of  the  pre -calculation  of  the  height  of  the  700-mb 
isobaric  surface  for  24  hours.  We  also  investigated  the  effect 
of  steric  spacing  upon  the  quality  of  the  prediction.  The  steric 
spacing  varied  only  with  the  calculation  of  the  laplacian  from  H, 
inasmuch  as  AH  is  the  most  sensitive  to  a  change  of  os  (deriva¬ 
tives  of  a  high  order).  The  spacings  Ss  =  250  and  As  =  500  km 
were  tested.  In  the  first  case  two  series  of  extreme  points  were 
lost,  and  consequently  H  was  determined  for  247  points;  and  in  the 
second  case,  three  series  of  points  were  lost,  and  H  was  deter¬ 
mined  for  187  points,  ' 

Finally,'  let  us  note  that  through  the  smoothing  of  the 
field  according  to  (10)  there  occurred  a  unique  correction  of  the 
values  of  H  in  the  external  zone  (excluding  the  external  contour) 
of  points,  for  which  H  is  not  calculated.  This  fact  can  be  con¬ 
sidered  as  a  certain  correction  of  boundary  conditions. 
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Hie  pre -calculation  of  the  height  of  the  TOO-nb  isobaric 
surface  was  carried  out  on  the  "Ural"  electronic  computer  in  the 
above -indicated  sequence.  The  program  was  placed  in  two  zones: 
the  program  was  entered  in  the  third  zone  and  the  original  data 
in  the  fourth.  A  block  diagram  of  the  program  is  presented  below 
/see  Figure  AppendixT" «  Selected  for  testing  the  efficiency  of  the 
prediction  system  was  the  period  14  through  18  November  1954,  which 
is  characterized  by  the  intensive  changeover  of  the  planetary 
height  frontal  zone  (PVFZ).  To  explain  the  potentialities  of  the 
baroti'opic  model  at  various  stages  of  the  transformation  of  the 
PVFZ  is  of  much  interest.  It  turned  out  that  fcr  the  given  period 
the  system  "cau^t"  best  the  development  (beginning  of  the  period) 
and  "caught"  worst  the  process  of  damping  (the  end  of  the  period) 
with  "erosion"  of  the  PVFZ  over  ETS,  For  various  time  (At)  and 
steric  (As)  spacings,  the  correlation  coefficient  between  the  cal¬ 
culated  and  actual  variability  of  the  heights  ranges  from  0,6  to 
0,75.  In  addition  to  this,  it  was  noted  that  the.  best  coincidence 
in  all  the  cases  which  were  analyzed  was  observed  for  As  «  ^00  km 
and  At  =  1  hour,  and  almost  the  very  same  thing  for  At  =  1,5 
hours.  For  the  steric  spacing  As  «  250  km,  the  correlation  co¬ 
efficient  for  all  time  spacings  was  at  a  minimumo  This  is  ap¬ 
parently  explained  by  the  fact  that  in  the  selection  of  such  a 
As,  the  defects  of  subjective  analysis  have  a  strong  effect  upon 
the  results  of  the  calculation.  With  the  selection  of  a  time 
spacing  At  >  1,5  hours,  there  appeared  spurious  waves  (first  of 
all  in  the  east),  which  reduced  the  quality  of  the  pre -calculated 
synoptic  situation. 

Examples  of  the  coincidence  of  the  calculated  (solid  lines ) 
and  the  actual  (broken  lines )  synoptic  position  are  presented  in 
Figure  la  (As  =  5^0  km.  At  =  1  hour),  lb  (As  «  5^0  hm.  At  = 

1,5  hours),  and  Ic  (As  =>  250  km.  At  =  1  hour).  For  the  original 
map  we  selected  AT -700  fc*r  18  hours,  l4  November  195^#  according 
to  which  we  calculated  the  iiap  AT-7OO  for  I8  hours,  I5  November 
1954  (i,e.,  the  prediction  for  24  hours),  Coiirparison  of  the  cal¬ 
culated  and  actual  maps  indicates  that  for  As  =  5 00  km  and  At  = 

1  hour  the  best  coincidence  occurs,  and  for  As  =  25O  km  and 
At  »  1  hOTTC*  the  worst  coincidence  occurs.  For  the  calculation 
of  AH  according  to  the  map  of  the  conventional  (subjective) 
analysis,  it  is  expedient  to  use  steric  spacing  of  the  order  of 
500  km.  For  the  determination  of  the  future  field  by  the  Euler 
formula  the  time  spacing  At  should  be  included  within  the  limits 
1  <  At  <2.  It  stands  to  reason  that  the  selection  of  the  time 
spacing  essentially  depends  upon  the  character  of  the  process 
which  is  occurring. 
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FIGURE  APPEMDK 


j  Entrance  of  the  program 
i  and  transfer  of  numbers 
jfrom  10  2 _ 


Block  diagram  of  the  program. 
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